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The properties of modes in a slab structure with gain in the center region 
and loss in the surrounding medium are investigated. The propagation 
constant and field distribution of the lowest-order modes are determined. 
The cutoff frequencies and propagation constants of the next-higher mode 
are given. Furthermore, the effect of a refractive index depression or in- 
crease in the center region is determined. The depression does not destroy 
the mode, as may be expected, but causes it to have a cutoff frequency. 
Comparison of these results to the experimental data shows that gain-induced 
modes play an important part in the lateral confinement in stripe- 
geometry GaAs lasers. 

I. INTRODUCTION 

Modes in cylindrical structures with refractive index boundaries are 
well known. Their existence and basic properties are generally visual- 
ized by superposition of plane waves reflected at the index boundaries. 
Much less is known about modes in structures with spatially non- 
uniform gain or attenuation. 12 However, it seems intuitively possible 
that some kind of mode should also exist in that case; indeed, Kogelnik 3 
showed that a cylindrical structure with a radial gain profile can sup- 
port a Gaussian beam of constant diameter even if there are no refrac- 
tive index differences present. Evidently, the nonuniform transparency 
of the medium counteracts the natural tendency of the beam to spread. 

In this paper we will consider a planar structure with stepwise dis- 
continuities of gain or absorption. This geometry is of considerable 
practical interest. As will be shown in this paper, the lateral confine- 
ment of the optical field in a stripe-geometry GaAs laser is due to the 
gain-loss interface at the edge of the stripe. Furthermore, it is easy to 
create nonuniform gain distributions in planar structures either by 
masking and optically pumping or by nonuniform injection current 
distribution. The modes induced by these gain distributions can be 
easily influenced from the outside by changing the pumping intensity 
or injection current. 
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II. GAIN-INDUCED MODES IN A THREE-LAYER STRUCTURE 

The physical configuration with which we will concern ourselves in 
this paper is shown in Fig. 1. All material parameters are assumed to 
be constant throughout the regions and differ only at the boundaries 
x = ± a. As is well known from the theory of refractive guiding, 4 the 
modes can be divided into four classes : TE even and odd in x direction 
and TM even and odd in x. The two lowest-order even modes have no 
cutoff frequency, whereas the odd modes can be guided only above a 
certain cutoff frequency. For our application, the lowest-order even 
modes are of greatest interest and we will deal with them first. 

Before going into details, let us explain the notation which we will 
use subsequently. 

(i) The complex relative dielectric constant e is split into refractive 
index n and extinction coefficient k according to 



e = (n + jk)\ 



(1) 



(ii) u determines the x dependence of the fields inside the center 
layer \x\ ^ a, i.e., they have the functional form 



os \ a ) 



sin 
cos 



(Hi) w gives the dependence of the fields in the cladding \x\ ^ a 

which follow the function e ~ wlxlal . u and w are related by 

W 2 + W 2 = ( €l _ C2 )(fc oa )2. (2) 

With these quantities we can derive the characteristic equations (the 



: 

a 

+a 


< 

z 


9 


/ 






—a 
a 





Fig. 1 — Cross section of slab waveguide. 
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field components are listed in Appendix A) for the even modes. 
TE Mode TM Mode 

w = u tan u — w = u tan u (3) 

They are formally the same as in the case of purely reactive modes but, 
since ei and e 2 are complex, the solutions will be complex. 

We determine under which conditions eqs. (3) have a solution. We 
restrict ourselves to cases where the fields cannot increase exponentially 
with increasing distance from the interface nor can there be any move- 
ment of wavefronts from infinity toward the guiding structure. This 
limits w to the first quadrant. It is then easy to derive a necessary but 
not sufficient condition from eqs. (2) and (3) : ei — e 2 has to stay in the 
first or second quadrant, i.e., Im (ei — e 2 ) > 0, which is equivalent to 

niki — n 2 k 2 > 0. (4) 

Let us interpret this inequality for some special cases : 

(i) ki = k 2 . Condition (4) reduces to n\ > n 2 which is the well- 
known requirement for refractive guiding. 

(w) nx = n 2 . The condition (4) now reads k x > k 2 , equivalent to 
saying that the center region should be more transparent than 
the sides, which agrees with physical intuition. 

(Hi) U\ = n 2 + An, where An is small compared to n 2 but not neces- 
sarily small compared to k. In this case (4) can be expressed by 

Ak , An ^ _ ,. N 

-j- H > 0, (4a) 

fci n 2 

where Ak = ki — k 2 . This inequality shows that there can be modes 
even if An or Ak is negative, as long as the other mechanism is strong 
enough to generate the mode. However, since (4) is only a necessary 
condition, this case requires further examination. 

The second problem we have to address ourselves to is that of the 
cutoff frequency. For refractive guidance the lowest-order even modes 
have no cutoff frequency. We will now establish the conditions under 
which gain-induced modes have no cutoff frequency either. For small 
values of (ei — e 2 )(fcoa) 2 , w and u will be small and the tangent in 
eq. (3) can be replaced by its argument. The solution of eqs. (2) 
and (3) is thus w — (ei — e 2 )(afc ) 2 for TE modes and w = («2/ci) 
X(«i — € 2 )(ak ) 2 for TM modes with ak very small. This result 
shows that the solution of (3) exists no matter how small (ak ) is, 
independent of the guide parameters, as long as ei — c 2 is in the 
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first quadrant. In other words, a gain-induced mode will have no cutoff 
frequency as long as there is a refractive index difference to keep ei — 62 
in the first quadrant. If there exists a refractive index depression, i.e., 
Re (ei — €2) < 0, the modes do have a cutoff frequency. 

We will now consider the case of gain-induced modes without a refrac- 
tive index change. We specialize the discussion to small extinction co- 
efficients in the order of 10~ 4 , whereas n is typically greater than one. 
This covers typical laser applications reasonabl3 r well. The character- 
istic equations for TE and TM modes [eqs. (3)] differ only by a factor 
of €i/t2 = 1 + j(2/n)(ki — k 2 ) which is sufficiently close to one to 
cause only a negligible perturbation. In the following, we will therefore 
neglect the difference between TE and TM modes. 

To simplify the discussion we define two new quantities : A normal- 
ized frequency is given by 

v = afcoVei — C2- (5) 

For gain-induced modes, v reduces to the form 



Sp {ki " hi) 



v = ak nJj2 ^ 1 * v , (6) 



which means that in this particular case the phase angle of v is inde- 
pendent of the material parameters. It is furthermore customary 5 to 
use a normalized propagation constant b 

o = or ( r- j = e 2 + (ei — €2)0. (7) 

tl — 62 \Kq / 

Since in our application ei — 62 is a small quantity, (/3,/fco) 2 is always 
equal to €2 plus a small perturbation, whereas the variation of 1 6 1 is in 
the order of one, thus alleviating some computational pro blems. The 
two parameters u and w are related to v and by u = Wl — 6 and 
w = Wo\ 

With these new quantities, the two characteristic eqs. (3) are trans- 
formed into 

V&" = Vn^ft tan (vjr=b). (8) 

It should be noted that the normalized propagation constant b is 
exclusively a function of v. It is, therefore, only necessary to solve the 
characteristic eq. (8) once to cover all material parameters and di- 
mensions, assuming, of course, the validity of the initial assumptions. 

The solution of the characteristic equation was done on the com- 
puter. Figure 2 shows the normalized propagation constant 6 for the 
loweet-order mode as a function of \v\. We can interpret b more easily 
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Fig. 2 — Z>(|y |) for lowest-order even mode (ni = n-i). 



if we rewrite the propagation constant fi z [eq. (7)] to show its real and 
imaginary parts. 



n - bi(ki - fc 2 ) + i(fc 2 (l - 6r) + fcl&r)- 



(9) 



The attenuation or gain of the mode as represented by the imaginary 
part of 0z/ko is solely dependent on b T , which varies between and 1. 
For small values of b r , the attenuation will be determined by the outer 
medium and only for higher values of b r will the gain in the center 
region be of significance. This situation corresponds very closely to the 
case of refractive guidance. For further reference we have included in 
Fig. 3 the field distribution with | v | as parameter. 

For any practical application, it is important to know up to which 
frequency or dimension the guide is single moded. We will therefore 
determine the cutoff frequency of the lowest-order odd mode and of the 
first-order even mode. We define the cutoff value of b as the one at 
which the radiation condition just ceases to be fulfilled. It turns out 
that this is the case when w is purely imaginary, i.e., only power is 
radiated away from the guiding structure, but the amplitude does not 
decrease with increasing \x\. Re (w) = is equivalent with the condi- 
tion Re (6) = [eq. (7)]. The solution of this problem has to be found 
on the computer. The results are 



fodd 



= 1.877 



= 2.759. 



For future reference we have included 6( | v \ ) for the first-order even and 
the lowest-order odd mode (Figs. 4 and 5). 
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Fig. 3 — Field intensity for lowest-order even mode. 



III. THE INFLUENCE OF ADDITIONAL REFRACTIVE INDEX DIFFERENCES 

In general, if there is a gain or loss difference between two regions, 
the refractive index will be different as well. Gain and refractive index 
are related by the Kramers-Kronig relations and thus the presence of 
gain can change the refractive index. Furthermore, in the case of the 
injection laser, the injected carriers change the refractive index due to 
the plasma resonance. In this particular case, the gain region can be ex- 
pected to have a lower refractive index than the side regions. It is thus 
necessary to explore the effects of refractive index increases and depres- 
sions in the center region in conjunction with gain-loss differences. 
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Fig. 4 — 6(|«|) for first-order even mode. 
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Fig. 5 — b{\v\) for lowest-order odd mode. 
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If there is an increase in refractive index in the center, we can expect 
the confinement of the energy to increase with increasing index differ- 
ence. However, in the case of a refractive index depression in the center, 
the two confining mechanisms will counteract each other. 

Let us consider the effect of increased refractive index in detail. In 
contrast to the case of gain-induced modes where the phase of v was 
t/4 [eq. (6)] independent of material parameters, it is now dependent 
on the complex difference ei — €2. This precludes a representation of 
the guiding parameters independently of the waveguide parameters as 
it was possible in the previous case. We are thus forced to determine 
a solution of the characteristic equation for each set of parameters. It 
is therefore of importance to develop approximations for the propaga- 
tion constant. We treat first the case, Ak/ An <K 1 . We regard the wave- 
guide as a perturbed refractive guide. The b(v) characteristic is plotted 
in Fig. 6. Obviously, both v and b are real. A small perturbation Ak has 
two effects: (i) the characteristic eq. (8) becomes complex and yields 
an imaginary part &,-, whereas the real part is to a first order unper- 
turbed; and (it) the propagation constant (9) is now approximated by 

p « n + Anbr -b<Ak + j(k 2 + b r Ak + b t An), (10) 

i.e., there appears a second imaginary component 6, An, which acts 
like a gain since An and bi are both positive. This effect is due to the 
improvement in guiding by the gain-loss mechanism. The calculations 
are straightforward and are listed in Appendix B. We note here only 
the result: 

bi = Ak (l - b )v W 
An 1 + yo V6o 

where the quantities inde xed wit h a zero denote the unperturbed state. 
In particular, Vo = ak n^2An/n. We notice that &,• contains a factor 
1 — 60, which decreases with b Q approaching unity. Inserting (11) into 
eq. (10) yields an expression for the loss (or gain) of the mode: 

*~4*'^ + <7tSt)]- 

We note that the function 

., . _ . (1 - 6 )t;oV6~o 
Kva) ~ bo + 1+..VS 
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Fig. 6 — b(v) for lowest-order even mode in a refractive slab guide and loss-weighing 
function f(v ) for lowest-order even mode. 



associated with Ak is only dependent on Vq and we have thus again 
achieved a representation valid for all waveguide parameters. Figure 
6 shows that f(vo) approaches unity much faster than 6o(vo); the loss 
is approaching the loss or gain in the center region indicating improved 
confinement due to the gain-loss difference. 

We will now consider the refractive index depression in the center 
region. From our discussion of cutoff frequencies, we know that modes 
in this case do have a cutoff frequency if they exist at all. It can be 
determined from the condition that w = v\lb is purely imaginary. 

In Fig. 7 real and imaginary parts of b are plotted as functions of 
s/X with the parameters listed in the figure caption. Since b r is larger 
than one, eq. (10) shows that the effective refractive index is smaller 
than that of either medium (An is negative). The phase velocity is 
therefore larger than that of a plane wave in either medium, as in 
metallic waveguides. Altogether, this is quite a different behavior from 
the ordinary dielectric waveguide. One would expect the index depres- 
sion to counteract the confinement and, if the depression is strong 
enough, to destroy it completely. As discussed before, this is only the 
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Fig. 7 — b(s/X) for lowest-order mode for a refractive index depression An = 10 -2 , 
n = 3.56, ki = 2.10-*, k 2 = 10"*. 



case below the cutoff frequency. To ensure that the solution of the 
characteristic equation is real and not a freak of the computer program, 
we have derived an approximate solution for Ak/An <5C 1 in Appendix 
C. The normalized propagation constant is given by 



6=1 + 



where 



sh'K i+ 'ife)}' 



v = - wn 



4 



2 1 Ant 



IV. GAIN-INDUCED MODES IN STRIPE-GEOMETRY GaAs LASERS 

The function of the stripe geometry has been viewed as selecting 
a filament and preventing others from forming. 6 We will now show that, 
in contrast to a laser with a wide area contact, the stripe geometry 
does provide a confining mechanism for the optical power and is not 
merely "selecting" the filament. 

In stripe-geometry lasers, the flow of carriers is confined in a stripe 
region parallel to the junction (laterally). This is done by proton bom- 



MODES IN PLANAR STRUCTURES 897 

barding the adjacent regions. 7 The proton bombardment and subse- 
quent annealing alter the free-carrier density but not the optical prop- 
erties of the material. 8 There are no deliberately built-in refractive 
index changes between the active and the passive regions. The experi- 
ments show that with stripe widths of ~12 ^m the lasers operate gen- 
erally in a single lateral mode; the origin of this confinement has not 
yet been explained. Figure 8a gives a linear scan of a typical nearfield 
distribution. However, at high current levels or wider stripe widths, 
numerous deviations from such a single-mode field distribution occur. 
Figure 8b gives an example of the nearfield in such a case. The field still 
fills essentially the whole stripe width but is highly nonuniform. Fre- 
quently the intensity maximum changes its position with changing 
current. It is clear that the simple case we have been analyzing so far 
cannot explain these effects. However, not enough knowledge is avail- 
able at this time about the numerous parameters involved. We there- 
fore try only to isolate the common properties of the majority of stripe- 
geometry lasers and explain them in terms of gain-induced modes. We 
had three major facts to consider: 

(i) For stripe widths in the order of 12 ^m and less, there is gen- 
erally a single mode for current levels close to threshold. 

(ii) At stripe widths above ~18 urn, the field distribution is very 
often nonuniform even at threshold. 

(Hi) The threshold current increases steeply with decreasing stripe 
width. 

Since we know that there must be a gain-loss difference between the 
stripe and the surrounding regions, we will apply the previous results 
on gain-induced modes and investigate if the experimental results can 
be explained by this effect or if additional mechanisms have to be in- 
voked. In doing so we must keep in mind that we can expect the 
agreement between experiment and this theory to be only qualitative 
since a number of effects are neglected. We mention the few most 
obvious : 

(i) The loss in the layers above and below the active region 
contributes to the overall loss. This contribution is dependent 
on stripe width. 
(ii) The gain depends on the field intensity and is therefore not 

constant. 
(Hi) The mirror loss may depend on the stripe width. 

The mirror reflectivity R of the lowest-order mode was determined to 
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Fig. 8 — Linear scan of nearfield distribution of stripe geometry lasers with stripe 
width of ~12 fim. One division =3.7 pm. Current in both cases is ~20 percent above 
threshold. 

be ~0.3 and hence the mirror loss, defined by a M = 1/L In l/R, to be 
^10/cm for L = 400 /xm sample length. 9 The attenuation of GaAs is 
~ 10/cm. We ask now how much gain would one need to reach thresh- 
old with an active region width of ~-12 faa. In the following, we will 
call gain g the excess gain over the intrinsic attenuation in the active 
region. To reach threshold, the effective gain of the lowest-order mode 
(gb r — a(l — b T ) from eq. (9)) must equal the mirror losses: 

gbr - a(l - 6 r ) = aM, (14) 
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Fig. 9 — Threshold gain as function of stripe width. Mirror loss ocm = 10/cm, 
n = 3.56, a = 10/cm, and refractive index depressions An = 10~ 3 , 10~ 4 . 



where a is the loss in the side regions. 6 r -as pointed out before-is ex- 
clusively a function of \v\ (Fig. 2). g\/^iv is plotted as a function of 
stripe width in Fig. 9. We can see that the gain value corresponding to 
~12 (Uin stripe width is quite realistic and we conclude that the gain- 
induced guiding must play at least some role in the stripe-geometry 
confinement. As a further piece of experimental evidence, it has been 
shown by Dyment 10 that the threshold current in a stripe-geometry 
laser increases strongly with decreasing stripe width. There is, unfor- 
tunately, some uncertainty about the relationship between current 
and gain. Generally, an exponential dependency is assumed (g ~ J", 
where q varies between 1.5 and 3). We have used q = 2 to insert 
Dyment's measured results into Fig. 9. We have furthermore intro- 
duced g(s/X) curves for refractive index depressions of 10 -4 and 10~ 3 
which correspond to carrier densities of 10 18 and 10 19 /cm, respectively. 
The agreement is satisfactory considering the accuracy of the measure- 
ments. It therefore seems reasonable to conclude that the gain-induced 
modes provide the dominant confining mechanism in stripe-geometry 
lasers. In view of this conclusion we will now derive a few properties of 
the gain-induced modes which should be useful for predicting some 



900 THE BELL SYSTEM TECHNICAL JOURNAL, JULY-AUGUST 1973 



properties of the stripe-geometry lasers. Wherever the ouput of the 
laser is to be coupled into an optical system, it is desirable to reduce the 
stripe width without increase in threshold current. Equation (14) was 
therefore evaluated for various loss values in the side region. Figure 10 
shows that-at least for zero refractive index difference—a has very little 
influence on the necessary gain. Evidently the increased loss in the 
higher a case is compensated for by improved guidance. 

A refractive index depression generates a gain-width relationship of 
the form [eq. (13)] 

, fl>0 i 1 /i c\ 

47m 2 (s/X) 3 /2Ari 



4 



It is quite evident that it is very difficult to compensate for the very 
strong s/X dependence by proper choice of An. It is, however, possible 
to reduce the stripe width with a refractive index increase, in which 
case the effective gain is given by [eq. (12)] 

a M + a(l - f(v)) 



/GO 



(16) 



where f(v) is plotted in Fig. 6. However, at present, the technological 
difficulties of this approach have not been solved. 
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Fig. 10 — Threahold gain as function of stripe widths for different attenuation a in 
the side regions, an = 12.5/cm. 
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V. CONCLUSION 



Modes caused by gain attenuation differences in a slab structure 
have been investigated. They can exist without refractive index differ- 
ences, in which case the lowest-order even modes do not have a cutoff 
frequency. The cutoff frequencies for the next-higher modes are given. 
The influence of additional refractive index increase or decrease in the 
center region is studied. One particularly interesting result is that, in 
the presence of gain, a refractive index depression does not remove the 
mode, but it causes it to have a cutoff frequency. 

An application of these results to the problem of lateral confinement 
in stripe-geometry lasers shows that gain-induced modes are mainly 
responsible for the optical confinement. Thus the stripe geometry does 
not just select a filament, but provides a confining mechanism for the 
optical power. (A similar conclusion was reached independently by F. 
R. Nash. 11 ) The sharp increase of the threshold current with reduced 
stripe width is a direct consequence of this gain-induced confinement. 
It therefore appears impossible to reduce the stripe width significantly 
below ~ 10 nm without a substantial sacrifice in the threshold current. 
A deliberately built-in refractive index increase in the center region 
should alleviate this situation. 

APPENDIX A 

Field Components in Slab Guide 

TE TM 

\x\ ^ a 

Ey = Auno cos u - e~ J 'P" H y = Aweoei cos u - e~ j&,z 

a a 

H x = —Apt cos u - e~ iP ' z E x = A@ z cos u - e~ j0 " 

a a 

H x = — jA - sin u - e~ i&tZ E z = jA - sin u - e~ i&tz 

J a a J a a 

\x\ ^ a 

E y = Aw/io cos we -W#.«+«(l*/«l-»] H y = iu6 e 2 cos we -'-''* / « z+u ' ( l* /a i- 1 >] 

H x = -Ap. cos ne-MfcHM.<i*/.|-i>] Ex = Ap z tl cos w -[**+.(i./.i-i)] 

H z = =FjA ™cosue-W'* + w"'\- 1 » E z = ±jA-^ 
a a 62 

Xcosue-^" +w ^ x '^- 1 ^ 
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The notations are explained in the list of notations at the end of the 
paper. The separation parameters are related to each other by 

OM 2 + u 2 = ei(froa) 2 , (0*a) 2 ~w 2 = € 2 (fc a) 2 , 

and u 2 + w 2 = (ei — e 2 )(fc a) 2 . 

APPENDIX B 

Approximate Solution of Characteristic Equation for Refractive Index 
Increase in Outer Region 

We assume that Ak/An <JC 1, i.e., the refractive index step is govern- 
ing the guiding properties. The normalized frequency v can then be 
expressed by 



{"fK' + 'E)- (17 > 



We treat the characteristic equation for TE modes first [eqs. (3)]. The 
regular solution for Ak = will yield real values for u and w. For 
small Ak/An, both will be complex with a small imaginary part. We 
list the real and imaginary part of the characteristic equation separately, 
neglecting second-order terms : 

u r sin u T = w r cos u r (18a) 

w,£(l + w r ) sin u T + u T cos u r ~\ = u>i cos u r . (18b) 

The real part (18a) is identical to the unperturbed equation and we as- 
sume, therefore, u r and w T to take the values of u and w, respectively, 
of the Ak = case. Combining (18a) and (18b) yields the following 
relation between w t - and «>,-: 

Ui\_Ur + Wr(l + Wr)] — WiU r = 0. (19) 

Now we have to relate Ui, Wi to 6,-, the quantity we really want to 
determine. 

w = Wb u = Wl - b. (20) 

With a small imaginary part of b we get 

■-*VK[i + >(£ + £)] (21a) 

.-^VT^Jn-^g-j^)]. (21b) 

These equations show that, if the real parts of w and w are unperturbed 
to the first order, the real part of b will be unperturbed also. If we intro- 
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duce eqs. (21) into the characteristic eq. (19), 6, can be determined 



= Afc (1 - b r )Vr^br 
An 1 + Vr^ r 



(22) 




The characteristic equation for the TM mode has an additional factor 
«i/e 2 . It is easy to show that the real part of b is altered by Ab from the 
TE value b 



(23) 



The imaginary part is to a first approximation the same as for the TE 
mode [eq. (22)]. 

APPENDIX C 

Approximate Solution of Characteristic Equation for Refractive Index 
Depression in Center Region 

Again we assume Ak/ \An\ <3C 1. An will be negative. The normalized 
frequency v is now expressed by 



.a $[En\(, . . Ak \ . /, . Afc \ 

V = 3 X n *V » { l + J 2|A^| ) = JVi { 1 ~ 3 2lAn~\ ) 



Let us assume that u is very close to ir/2, 

u = tt/2 + 8 + jf, 
where 5 and \p are small quantities. Since 



(24) 



(25) 



w = 



we use 



w = v 



-mi 



(26) 



where we note that u/v is small compared to one. The characteristic 
equation to a first-order approximation takes the form 



(| + * + #) +**(« + #) =°- 



The solution is 



• --; ! 



21 + »? 



+ = l Vi 



2 1 + v\ 



(27) 



(28) 
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It is now easy to calculate b, which is related to u by 



The result is 




(30) 



keeping in mind that v r /v, <5C 1 and v \ » 1 . 



NOTATIONS 




a = 


47rfc 
X 




a 






b 






s = 


2a 




An 


= Wi — 


n 2 


Ak 


= fcx- 


fc 2 


to 






«1» « 


2 




9 






/Co = 


: O)V/X0€0 


fcl,2 






X 






Mo 






Wl,2 







attenuation constant 

halfwidth of slab 

normahzed propagation constant [[defined in eq. (7)2 

width of slab 



dielectric constant of vacuum 

relative dielectric constants (subscript 1 refers to 

the center and 2 to the outer region) 

gain in center region 

free-space wave number 

extinction coefficient (subscript 1 refers to the center 

and 2 to the outer region) 

free-space wavelength 

permeability of free space 

refractive index (subscript 1 refers to the center and 

2 to the outer region) 

angular frequency 

u separation parameter 

v = afcoVei — € 2 normahzed frequency 
w separation parameter 
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